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We examine static perfect fluid spheres in the presence of a cosmological constant. Due to the 
cosmological constant, new classes of exact matter solutions are found. One class of solutions 
requires the Nariai metric in the vacuum region. Another class generalizes the Einstein static 
universe such that neither its energy density nor its pressure is constant throughout the spacetime. 
Using analytical techniques we derive conditions depending on the equation of state to locate the 
vanishing pressure surface. This surface can in general be located in regions, where going outwards, 
the area of the spheres associated with the group of spherical symmetry is decreasing. We use 
numerical methods to integrate the field equations for realistic equations of state and find consistent 
results. 
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^ ^ I. INTRODUCTION 

^ ' Static and spherically symmetric perfect fluid solutions have always been a rich source of investigation in classical 

5-H , general relativity ever since the pioneering work of Schwarzschild in 1916. He solved the field equations for the interior 

I ^^ ' region by assuming a perfect fluid of constant energy density, and also for the outside vacuum region, the famous 

Schwarzschild solutions. 

C\| ■ The interior solution has the geometry of a three-sphere, which was noticed as early as 1919 by Weyl p|. The 

^ ' pressure of the Schwarzschild interior solution always vanishes before the equator of the three-sphere, therefore, this 

^^ . geometrical picture is not of great importance. This situation chances significantly in the presence of the cosmological 

constant [3, y|- In some special cases the interior metric with cosmological constant was studied earlier, see e.g. [J, Q- 

However, the complete analysis was only completed Q, y] 85 years after Weyl noted the interesting geometrical 

structure of these solutions: With A the pressure can vanish exactly at the equator of the three-sphere in which case 

one has to join on the Nariai metric p, |7[ as the exterior vacuum metric. We note that although the spatial geometry 

of the interior Schwarzschild is a three sphere, the four-metric is not homogeneous. There is a center of symmetry in 

i"-^ ' the fluid region, surrounded with concentric spheres defined by the orbits of the group generating spherical symmetry. 

• • . The spacetime metric associate an induced metric and thereby an area to each group orbit. Going outwards from 

the center this area is increasing until one reaches the equator, the sphere with the maximal area. It is furthermore 

possible that the pressure vanishes in a region where the area of the group orbits is decreasing. 

^ • Lastly, the matter can occupy the whole three-sphere having two regular centers. This generalizes the Einstein 

5^ \ static universe 0, y, la] ' see @ for early results in that direction. The Einstein universe also emerges as a special 

case (vanishing expansion rate and vanishing vorticity) when considering homogeneous shear-free perfect fluids [lO[. 

In recent years it has been generalized in various different theories, like brane world models |ll| . Einstein-Cartan 

theory [l2| , in modified gravity theories [iJl or Loop Quantum Gravity [ij] . 

In this paper we are analyzing systematically the effects of a positive cosmological constant on perfect fluid spheres. 
We generalize known exact solutions of the field equations with cosmological term and discuss their new properties. 
The principal result of the analysis is the increase of the radial sizes of matter spheres. This naturally lead to questions 
regarding the physical picture applicable to these solutions. Analytical and numerical techniques are used to obtain 
a consistent picture of the underlying physics. 

A positive cosmological constant can be regarded as an external force pulling matter apart. Therefore, a 'large' 
positive cosmological constant can increase the radius of a known perfect fluid solution such that it occupies more than 
just 'half of the three-space. The effects of the actual cosmological constant are very small and therefore one may 
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ask why the investigation of solutions with 'large' values of the cosmological term is beneficial. This can be answered 
from two points of view. From a mathematical point of view we deepen our understanding about exact solutions of 
Einstein's field equations and the influence of the additional parameter A. However, also from a physical point of 
view this study can be justified easily. In the bag model of hadrons [la, [l^, \1^ the bag is stabilized by a term of the 
form Bgab, which has the same form as a cosmological constant, though its numerical value is considerably bigger 
Be ~ 8.91MeV, whereas we have A^'"* « 1.78 x 10~^MeV. Hence, the key physical motivation is the possibility that 
effects within stellar models can be effectively described by a term that looks like a cosmological constant, which can 
indeed have large effects. Also in the context of Boson stars jlSl ll9| the matter energy-momentum tensor contains 
a part proportional to the metric which can be read as an effective cosmological term. Within the context of Loop 
Quantum Gravity it has recently been shown [20[ that quantum gravity effects can be effectively of the form of a 
cosmological constant. 

Apart from an equation of state relating the density p to the pressure p, a spherically symmetric perfect fluid 
has to satisfy only one condition, the pressure isotropy condition, requiring the equality of the radial and angular 
directional pressures. Since this condition does not involve the cosmological constant A, any solution with A = 
having an equation of state f{p,p) = can also be interpreted as a A 7^ cosmological solution with equation of 
state f{p + A,p — A) = 0. It is one of the main purposes of the present paper to re-investigate known perfect fluid 
solutions with cosmological constant and to analyze whether these solutions for some special cosmological constant 
require the Nariai metric as the exterior vacuum spacetime. We also try to construct solutions with two regular 
centers, i.e. constructing more general Einstein static universes, having neither constant energy density nor constant 
pressure. 

Since the static and spherically symmetric field equations with more realistic equations of state in general cannot be 
integrated analytically, we also use numerical methods to study perfect fluid spheres. In particular we will show that 
the effect of a 'large' cosmological constant on polytropic perfect fluids is such that the matter is pulled sufficiently 
apart so that it occupies much of the three-space. The same result is also found by considering the stiff matter equation 
of state and also the Hagedorn equation of state. However, realistic equations of state seem not to allow the presence 
of a second regular center. It is important to distinguish between coordinate and physical effects. For instance, one 
cannot numerically integrate the constant density solutions for large cosmological constants if the original radius is 
used as a variable. The code breaks down at the equator of the three-sphere. This is a pure coordinate effect, since 
the coordinate system does not cover the whole spacetime. In order to distinguish coordinate and physical effects, 
also the Ricmann curvature tensor is considered. 

Various other effects of the cosmological constant have been studied in the past, like the dynamical instability of 
perfect fluid spheres [2l|, [24I , possibly detectable effects of the cosmological constant within our solar system [2^, ^M, 
\2dl . [26| , and effects within astrophysical structures [23, [2^, [23, [30, [3l| . Recently, the bending of light with A has been 
discussed in [33, [33| . 

The paper is outlined in the following manner: In Section [Hj we analytically study the effect of the cosmological 
constant on the Whittakcr and Tolman IV solution, and also discuss the matching of the matter and the vacuum 
solution. In Section [TTTl we derive conditions on the equation of state and the cosmological constant to characterize the 
different possible solutions. In Section [TVl wc numerically integrate the field equations for realistic equations of state 
and find results consistent with our analytical results. We summarize and conclude our work in the final Section |V| 

II. THE WHITTAKER AND TOLMAN SOLUTIONS 

In the present section the Whittaker solution [3J] and the Tolman IV solution [33 are recalled. In both cases we 
firstly introduce a third angular coordinate a, so that the coordinate system covers the complete three-space and not 
just 'half of it like the usual radial coordinate r. Secondly, wc introduce the cosmological constant in the solutions 
and analyze its effect (the 'external' force due to A) on those solutions. We also explicitly show, how to join the 
interior and the exterior solution through the vanishing pressure surface. 

A. The Whittaker solution 

The Whittaker solution is characterized by the relation p + 'ip = po between the energy density and the pressure, 
where po is a positive constant. Similarly to the constant density case this condition allows one to write down the 
solution to Einstein's field equations in terms of elementary functions. The metric of the Whittaker solution in the 



original Schwarzschild coordinates reads [3j 



ds^ = -b 



B 



1 + B v 1 — a^r^ arcsin(ar) 



dt^ 



B 



(1 + BMl - a^r^) (1 - aV^)^/^ arcsin(ar) 

ar 



dr'^+r^dfl^, (1) 



where a, b and B are constants. The special importance of the Whittaker solution lies in the fact that it is the 
non-rotating static limit of the Wahlquist solution [36| . the most important rotating perfect fluid exact solution. The 
parameter k of the Wahlquist solution is related to the parameter B of the Whittaker metric hy B ~ 1/k^. Since the 
change of the parameter n corresponds merely to a rescaling of the coordinate t, we set 6=1. After introducing a 
new radial coordinate, the third angle a, by r = (1/a) sin a, metric ([1} simplifies to 



ds' 



-fdt^ 



1 /da^ 
T 



sin a dri" 



f = 1 + B{1 -a cot a). 



(2) 
(3) 



Although the introduced radial coordinate a is very similar to the third angle of the ellipsoid used for the interior 
Schwarzschild solution, the spatial metric is not ellipsoidal in this case because of the non constant nature of the metric 
component gaa- On the other hand, it should be emphasized that the r = const, hypersurfaces are round two-spheres. 
Therefore, the topology of the three-space essentially depends on the function gaa- In the above mentioned constant 
density case, the three-space is in fact a three-sphere. When we discuss next the modified Tolman IV solution the 
three-space will be ellipsoidal. The introduction of the new radial coordinate a is important in all cases where there 
is a group orbit with maximum area, since in these cases the usual radial coordinate only covers the region up the 
this maximum orbit. Henceforth we will refer to the coordinate a as the third angle. 

It is possible to express the constant in the equation of state p + 3p = po in terms of the constants in the metric 
and the cosmological constant by 



Po 



a^S + A 

An 



Positivity of the pressure and density implies po > 0; which we require from now on, by assuming 

A > -a^B. 



(4) 



(5) 



If there is a spherical surface where the solution is matched to a Schwarzschild-de Sitter (or Schwarzschild anti-de 
Sitter) exterior region, then p must go to zero at the surface, and the fluid density becomes po there. 
The pressure and energy density of the Whittaker solution are given by 

(6) 

(7) 

The central pressure and central energy density can be derived from equations ([6]) and ([7]) by noting that liniQ^o ./ = 1 
and read 
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(8) 



Requiring both to be positive we get ^po < a^ < AirpQ. Since at the center 



df 



and 



d'/ _ 2B 



-j^ = ^, the pressure is 

maximal at the center if and only if B > 0. Since in the B < case there is no zero pressure surface, we assume 
i? > in the present discussion. 

The equator, where the area of the spheres of symmetry is maximal, is located at a = n/2. There f = I + B and 
for the equatorial pressure peq and density peq we get 



Peg = -^{^~a^), 



Peq 



-_ia^^A) + PQ. 



(9) 



Let us now choose the cosmological constant such that the pressure vanishes at the equator, p^q = in Q, which 
yields 

A = a^ =: An. (10) 

Hence for solutions which satisfy the relation A < A^r, the pressure vanishes before the equator, where the group 
orbits are still increasing. 

One should check whether the choice A = An is compatible with the positivity of pressure and energy density 
at the center, and also with the positivity of energy density at the equator (otherwise these solutions would not be 
physical). Positivity of energy density at the surface is ensured by ([5]), which now takes the form B > —1. Using 
A = a^ in (g]) and ^ we get 

p, = ^B> 0, (11) 

p, = ^(2-i?)>0. (12) 

All the three conditions are satisfied if < i? < 2. 

For larger cosmological constants, i.e. for A > Ajv, the pressure vanishes after the equator of the ellipsoid, where the 
group orbits are decreasing. Since the pressure function p —^ — oo as a —^ tt there always exists a zero pressure surface 
at which one joins on the Schwarzschild-de Sitter metric as an exterior vacuum spacetime. Therefore the Whittaker 
solution cannot have a second regular center. Since the derivative of the function / is positive for < a < tt in the 
case B > 0, the function / remains positive, consequently the coordinate system and the metric remains regular in 
the whole fluid region. 

As already outlined in the introduction we will explicitly show that the choice A = Ajv necessitates the Nariai 
metric as the exterior spacetime. Since the pressure in this case vanishes at the maximum of the area of the group 
orbits the exterior spacetime must have constant area spheres of symmetry, which excludes the Schwarzschild-de Sitter 
or anti-de Sitter spacetimes. The only other static and spherically symmetric vacuum spacetime with cosmological 
constant is indeed the Nariai spacetime, and its group orbits have constant area. 

B. Matching procedure 

Here we review the necessary conditions for matching a spherically symmetric static perfect fluid solution to an 
exterior vacuum region. We write the metric in both regions in the form 

ds^ = -e^^df + \dr^ + R^{de^ + sin^ 6'd<?!)2), (13) 

y 



where v, y and R are functions of the coordinate r. A metric written in this form have a center where i? = 0, and this 
center is regular, i.e. free of conical singularities, if the area of small spheres is proportional to their radius square, 
with the appropriate proportionality factor, ^ = ±i. The regularity of the four-metric also requires a finite value 

for V and ^ = at the central point. 

We assume that the matching is performed along the hypersurface described by r = Ts- The Darmois-Israel 
matching conditions [33, [33| essentially state that the induced metric and the extrinsic curvature have to agree on the 
hypersurfaces used for joining the two solutions. The outward pointing normal vector to the symmetry surfaces has the 
components n° = (0, y, 0, 0) . The induced metric hab can be expressed using the spacetime metric as hab = gab — naTib, 
while the extrinsic curvature can be calculated as Kab = h^V cnb^ where Vc denotes the covariant derivative associated 
to the spacetime metric gab- Using the coordinate system {t,9,(j)) on the matching surface, the induced metric hab 
has the form 

e^-^ \ 

hab=\ R^ , (14) 

i?2 sin^ 9 J 

while the components of the extrinsic curvature are 

-ye^'^f- 

Kab=\ ' yR^ I . (15) 

yR^ sin^ 



From p4|) it is apparent that the induced metric agrees if and only if the value of R and i/ agrees on the matching 
surfaces. Then the extrinsic curvatures are matched appropriately if and only if y^ and y^ agrees. The agreement 
of R has the obvious physical meaning of equal area matching spheres, while the equality of v can always be ensured 
by appropriate rescaling the time coordinate in the interior fluid domain. If we use Gauss coordinates in both domain 
then y = 1 and the matching of the extrinsic curvature is equivalent to the continuity of the first derivative of v and 
R. However, in general, it is also possible to give coordinate system invariant meaning to these conditions. 
The invariant mass function in spherically symmetric cosmological spacetimes can be defined as 

m=|(l-5"'i?.ai?,6)-|i?'. (16) 



For vanishing cosmological constant this gives back the usual mass definition given in [39j . For the metric form ([Te 
we get 



R 

"2 



2 I diV 
dr 



^R'- (17) 



From this we can see that if m agrees on the two matching surfaces of equal area then y^ must also be the same. 
The other invariant quantity is the pressure at the matching surface, which for the metric (|13p takes the form 



y"^ dR fdv 1 dR\ 

It is apparent that if R and y^ both agree on the matching surfaces and ^ 7^ then y^ will be the same if and 
only if the pressures are the same. Consequently, if ^ 7^ 0, the matching of two static perfect fluid solutions can be 
done at two chosen spherical surfaces if and only if the surfaces has the same area, the mass function has the same 
value, and the pressures agrees as well. Obviously, if the exterior domain is a vacuum, then the fluid pressure at the 
surface must vanish. It is interesting that in case of a Nariai exterior -tt = and the p ~ Q condition is not enough 
to ensure the continuity oi y^. 



The quantity y4^ is closely related to the acceleration of static non-rotating observers staying at constant radius 



r. 



In the coordinate system x" = {t, r, 6, (f)) used in (fT3|) these observers have the velocity vector v°- — {e '^, 0, 0, 0). The 
only non- vanishing component of their acceleration vector a° = u^Vbu" is a'' = y^^. The norm of the acceleration 
is \a\ = ^/oFa^ ~ ^ |^|- This shows that apart from a possible signature change the continuity of the magnitude of 
the acceleration implies the continuity of y4^ in the matching condition. 



C. Joining interior and exterior solution 

For cosmological constants satisfying A < A^r the area of the group orbits at the p = Q surface is increasing and we 
join the Schwarzschild-de Sitter (or Schwarzschild anti-de Sitter for A < 0) metric on as the exterior vacuum spacetime. 
Since the cosmological constant is fixed by the specific solution it remains to choose the mass appropriately. In the 
Schwarzschfld area coordinate R the mass is defined by M = /^ ° 4:TrR^p{R)dR, where Rs is defined by p{Rs) = 0. 
By using Gauss coordinates relative to the r = const, hypersurfaces the metric is C^ at the boundary. If the energy 
density is non-vanishing at the boundary this cannot be improved. After placing one object in the Schwarzschild-de 
Sitter spacetime, it still contains an infinite series of singularities. However, by placing a second object appropriately 
in that spacetime, it is possible to construct a singularity- free spacetime, see Fig.[T] This possibility has been discussed 
earlier in greater detail in [3, Q • 

For A ~ An wc explicitely show the matching of the interior perfect fluid spacetime with the exterior Nariai 
spacetime. We follow the generic discussion of matching two static and spherically symmetric regions presented in 
the previous subsection. We read off and compare at the matching surface the corresponding functions v, y and R in 
the general form of the line element (fTH]) . 

Recall that the static form of the Nariai metric is given by [3] 

ds^ = -cos^xdt'^ + j(dx'^+dn'^). (19) 

We note that this form of the metric is not homogeneous, since the static observers described by constant (x, 9, 4>) 
are not equivalent. The magnitude of their acceleration is \a\ = vAtan^. Somewhat surprisingly, this acceleration 
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FIG. 1: Penrose-Carter diagram with two stellar objects separated by a Schwarzschild-de Sitter vacuum domain. The radii R 
of the stellar objects is between the radii of the black hole and cosmological horizons. Since the group orbits are increasing up 
to R the vacuum part contains the cosmological event horizon r++. 



is towards the x = surface, since the non-vanishing component of their acceleration, given by a'' 
negative for x > 0. 

The metric functions of the Nariai spacctime in the coordinate system (|13p are given by 



e^''" = cos2 X, y% = A, 

The corresponding functions in the Whittaker fluid region are 






-Atanx, is 



(20) 



r.'^'^W 



c'f, 



2 2 p 



R 



sin^ a 



w 



where / is given by ^. The agreement of the induced metric, i.e. the matching of R and v imphes 



1 

A 



sin a 



cos x^c-f. 



(21) 



(22) 



The condition yN^^§^ = yw '^^^ implies a ~ 7r/2 for the matching surface in the fluid region. From this it follows 
that A — a^ , which is just the condition of vanishing pressure at the equator of the Whittaker solution. The remaining 
matching condition yN^j^ ~ Vw^^^ yields 



tanx 



ttB 



_7r/ _47rpo\ / A 
4VTTB 4 ^ A y V Airpo ' 



(23) 



So if Gauss coordinates are used, i.e. y = 1, we explicitly showed the matching of the interior and the exterior metric 
to be of degree C^ (i/, R, v' and R! all agree on the zero pressure surface). This differentiability condition cannot be 
improved when the equation of state of the fluid constrains the energy density at the boundary to be positive. In this 
case the energy-momentum tensor jumps at the boundary and the metric is at most C^ . Figure [2] shows the Penrose 
diagram of this spacetime. 

Lastly, for cosmological constants satisfying A > Ajv the group orbits at the p = surfaces are decreasing and 
we join the Schwarzschild-de Sitter metric on as the exterior vacuum spacetime. Since the cosmological constant is 
fixed by the specific solution it remains to choose the mass appropriately, see the discussion above. By using Gauss 
coordinates the metric is C^ at the boundary. If the energy density is non-vanishing at the boundary this cannot 
be improved. It is important to note that the vacuum part of this spacetime contains the singularity at the origin 
r = and that the matter occupies the 'outer' region of the spacetime, see Penrose diagram[3l where a second object 
was inserted in the spacetime to remove the infinite sequence of singularities present in the Schwarzschild-de Sitter 
diagram. 



D. Cosmological Tolman solutions 



In principle, the above discussion can now be repeated for all the 127 candidate solutions presented in [40] (of 
which only 60% are isotropic and regular at the center). For all those one could check whether the inclusion of 




FIG. 2: Penrose-Carter diagram with two stellar objects having radii 7? which require the Nariai spacetime to be the vacuum 
part of the global solution. The solid and dashed lines represent the future and past event horizon, respectively. 
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FIG. 3: Penrose-Carter diagram with two stellar objects connected by a Schwarzschild-de Sitter vacuum domain. Since the 
group orbits are decreasing at the matching surface the vacuum part contains the black hole horizon and the r = singularity. 



the cosmological constant can pull the matter up to or beyond the equator of the corresponding ellipsoid. We will, 
however, only repeat this analysis for the Tolman IV solution. This particular choice is motivated by the simplicity 
of the solution which also allows to analytically express its equation of state. 
We now take a fresh look at the Tolman IV solution [Sal. Its metric reads 



ds^ = -bHi 



'■/A^)df 



1 + 2rVA2 



-dr^ 



(l-r2/i?2)(l + r7A2) 
where A and R are positive constants. Introducing the third angle by r = Rsina yields 

■1 + 2{R^/A^)sin^a ^ . 



dn^ 



ds' = -B^l + (i?VA^) sm'a)dt' + i?^ 
Pressure and energy density respectively are given by 

i?2 _ 3i?2 gij^2^ _ ^2 



1 + (i?2/A2) sin^a 



-da' 



'adn' 



8irp{a) 
8np{a) 



i?2(A2+2i?2sin2a) 

i?2 + 3i?2 gjj^2^ _^ 3^2 



A, 



2^2 cos2a 



i?2 ( A2 + 2i?2 sin'a) ( A2 + 2i?2 sin^a)2 



A. 



(24) 

(25) 

(26) 
(27) 



Since all quantities depend on a only through sin a and cos2 a, the solution is symmetric to the equator a = 7r/2. 
Eliminating the variable a from Eqs. (p6)) and (|27|) leads to the following equation of state 



P = C0 + C1P + C2P , 
where the three constants c^ are given by 

(3 - 2Ai?2)(i?2 + ^2(2 _ Ai?2)) 



Co = 



47ri?2(^2 _^ 2i?2) 



Cl 



2i?2 + A2(13-8Ai?2) 
A2 + 2i?2 ' 



C2 



327rA2i?2 
^2 + 2i?2 • 



(28) 



(29) 



As in the previous discussion let us now compute the value of the pressure and density at the first center (a = 0), 
at the equator of the ellipsoid (a = 7i'/2) and at the second possible center (a = tt), which yields 

8^p(0) = -1 - -1 + A, SnpiO) = A + ^ _ A, (30) 

8^p(V2) = -;^+A, 87:pin/2) = ^^^^^^-^^^-A, (31) 

^M^) = i - ;^ + ^' 87rp(^) = -1 + A _ A. (32) 

Similarly to the Whittaker case there exists a cosmological constant such that the pressure vanishes at the equator 

A=-A=:Aw, (33) 

where the energy density is positive. In this case one has to match the Nariai metric. For smaller cosmological 
constants, i.e. A < An, the pressure vanishes before the equator and one has to join on the Schwarzschild-de Sitter (or 
Schwarzschild anti-de Sitter) metric as the exterior vacuum spacetime. However, the pressure cannot vanish after the 
equator, as can be seen from the mirror symmetry to the equator of the solution. Therefore, solutions with A > Ajv 
have two centers. One easily verifies that both centers are regular by checking that the derivatives of pressure and 
energy density vanish at both centers. Hence, as a side result we already found a new generalization of the Einstein 
static universe. 

The Einstein universe is characterized by its constant energy energy and constant pressure (originally Einstein 
assumed a pressure-less universe) throughout the three-sphere. By a generalization of the Einstein static universe we 
mean a globally regular solution of the field equations with cosmological constant where the spatial part of the metric 
is a closed three-space and where either the energy density or the pressure, or both are varying. 

It is expected that other known perfect fluid solutions, e.g. those given in Ref. [40[, will show similar properties. 
Therefore, following the above procedure, we can explicitely show the existence of a wide class of generalized Einstein 
static universes and an even wider class of static and spherically symmetric perfect fluid solutions, for which the 
pressure vanishes in regions where the group orbits are decreasing. 

III. ANALYTIC CONSIDERATIONS 

In the previous sections we analyzed perfect fluid solutions which may extend through the equator of the ellipsoid 
that describes the global geometry of the spatial hypersurfaccs. The present section supplements the explicit and 
later numerical results by presenting some general statements. It should also be mentioned that the existence and 
uniqueness of perfect fluid solutions was proved in [4l| . The restrictions on the equations of state could be weakened 
in Refs. [43, |4^. The existence and uniqueness proof of [41| could be extended to include cosmological constants 
satisfying A < 47rp(p = 0) in @, |4J|. Let us now consider the static and spherically symmetric line clement in Gauss 
coordinates relative to the r = const, hypersurfaces 

ds^ = -e^^'-'^Ut^ + dr^ + R'^{r)dn^, (34) 

the resulting field equations Gab + Agat ~ STrTah are given by 

1 - i?'2 - 2RB!' 



i?2 
i?'2 - 1 + 2RB!v' 

v'B! + i?(z/'2 + v") + i?" 



- A = Sttp, (35) 

+ A = Sttp, (36) 

+ A = Sttp, (37) 



R 

that imply the conservation of the energy-momentum tensor 

p' + i^'(p + p) = 0, (38) 

where the prime denotes differentiation with respect to r. The maximum of the area of the group orbits is located at 
the maximum of the function R{r), which means R'{r„i) = 0, where r^ is the location of the maximum. We henceforth 
assume that such a maximum exists (To be precise we only assume that the function R{r) has an extremum and it 



will turn out that this cxtremum is a maximum if wc require the energy density to be positive at r^ . For r = r„j the 
mass definition (fT7|) implies 1 — 2m{rm) / R{rm) — ^/'!>R(j'mY — 0)- 
Note that mass (|17|) and energy density p are related by 

dm m' ^ , , 

^ = ^=4.pi?^. (39) 

Eliminating the function v' from the first two field equations (P5)) - ([55|) and the conservation equation ([55)1 yields 
the Tolman-Oppenheimcr-Volkoff (TOV) [3^, |43 equation 

dp ^ (p + p)(47rp + m/i?3-A/3) 

JR = -^ l-^-A^3 ' (40) 

where we used that dp/dr = {dp/ dR)R' . At the maximum r„i the TOV equation is ill-defined since the denominator 
tends to zero. However, this is not a physical singularity, as can easily be seen by considering the derivative of the 
second field equation (|M)) evaluated at i?' = which reads 

8V(r™)-2^'(r™)^^^, (41) 

R\rm) 

and moreover all Riemann tensor components (jAip arc well defined at r = r^- Furthermore one can the express 
energy density plus the pressure function in terms of the Riemann tensor, 

4^(p + p)= R-\.e - R^'et = R''^r^ - R'^'<pt, (42) 

so that this sum is well defined if the spacetime is non-singular. Let us furthermore evaluate the field equations at 
r„i which yields 

^-^^(-)^"(-l-A^8..(.„), (43) 



RHrm) 

R{rm){^'{rm? + i^"{r^)) + R"{rm) 
R{rm) 



R^rm) 

+ A = 8^p(r„), (44) 



+ A = 8^p(r„). (45) 



Next, from equation (|44p we find that the pressure is positive at the equator if the cosmological constant is large 
enough, this means if 

For the special case A = Ajv the pressure vanishes at the equator, i.e. the maximum of the area of the group orbits. 
Putting this particular value of the cosmological constant into the first field equations yields 

from which we conclude that R"{rm) < to have a physically meaningful perfect fluid solution. Indeed, this condition 
simply states that the equator is a local maximum of the group orbits' area. This fact was not assumed explicitely 
and is a direct consequence of assuming physical solutions: positivity of the energy density. 

The above analysis clarifies under which conditions static and spherically symmetric perfect fluid ellipsoids may 
extend through the equator and have a zero pressure surface in regions where the area of the group orbits is decreasing. 
However, condition (j46p depends on the function R(r) and therefore essentially depends on the solution. Given an 
equation of state, a central pressure and a cosmological constant, such that the pressure (and by the equation of state 
the energy density) is decreasing near the center, the above considerations are not sufficient to decide whether the 
pressure vanishes before the equator or not. Therefore one should find a condition depending solely on the equation 
of state and on the initial conditions (pressure at the center and the cosmological constant) so that one controls the 
location of the zero pressure surface. 
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The TOV equation (|40|) seems to be ill-defined at r^, since its denominator vanishes. On the other hand, the 
field equations imply that the spacetime is regular where R'^r^n) = 0. Therefore we can conclude that the limit 
lini,.^r„, p'if) rnust exist, so that we write linir_,r™ p'{f) = <i- Existence of p'(rm) can be put back into (|40|1 and yields 

^ p(r-„.)+p(r,„) (4^p(r)i?(r)3 + m{r) - A/3R{r)^) 

" Rir„,y r™„ R'{r) ' ^ ^ 

Since limr^r„ R'{f) ~ 0, the numerator also must vanish as r ^ rm, 

lim {4np{r) + m{r)/R{rf - A/3) = 4^p(r„) + m[r^)l R{rmf - A/3 = 0, (49) 

r — >-r^ 

which can be written more conveniently (for the present purpose) 

47rp(r„,) = j^^ (^^R{rraf - m{r„,)^ . (50) 

Before exploiting the latter equation ([5(I)l . we note that it is easy to show that the pressure in the TOV equation (PO]) 
is decreasing near the center if 

A < 47rp(pc) + 127rpc, (51) 

a condition that only depends on the initial values and the equation of state. 

According to equation (j50p the signature of the pressure at the equator is determined by the signature of the 
quantity 

A f^" 

-i?i - m{R„,) = [A - ATip{R)]R^dR, (52) 

•J Jo 

where we used the definition of mass, sec (|39p . Since pressure and by a monotonic equation of state also the energy 
density are decreasing functions, p{p = pc) > p{rm), and a sufficient condition to satisfy the inequality p{rm) > is 

4ttp{p = Pc) < A. (53) 

Hence, if we the cosmological term is large enough, compared to the central density, then the pressure does not vanish 
before the equator of the interior spacetime. 

Next we find a necessary condition for a positive pressure at the equator. Let us assume p{rm) > 0. Then the 
integral of A — AnplR) is non negative. Since going outwards p and p are monotonically decreasing it follows that 
A — 47rp(i?„i) > 0. But because the zero pressure surface is after the equator, by the monotonicity condition we have 
p{Rm) > p{p = 0); and consequently A > 47rp(p = 0). On the other hand, this means that if 

A < 47rp(p = 0) (54) 

then necessarily p(rm) < 0. This condition is in agreement with previous results, see e.g. 0,3- Hence, if the given 
equation of state and the cosmological constant satisfy the condition (|M1) . then the pressure vanishes before the 
equator of the ellipsoid. 

Next, let us assume that the pressure vanishes at the maximum of the area of the group orbits, so that the equator 
is also the zero pressure surface. In this case, as we already showed in the previous sections, one has to join on the 
Nariai metric. Putting p{rm) = p{rb) = into ((501) leads to 

m{n) = M = -Rl [ \^p{R)dR=-Rl (55) 

which relates the mass of the solution to the cosmological constant. Unfortunately this condition cannot be written 
in a form such that the equation of state suffices to choose the initial values for such a solution. 
These three observations can be summarized as follows 

A < Attp{p — 0) pressure vanishes before the equator, 

A'Kp{p = 0) < A < 4ttp{p = Pc) no analytical control, 

4:'Kp{p — Pc) < A pressure can vanish only after the equator. 
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which only depend on the equation of state and the initial conditions. However, these conditions only yield quite 
general conclusions. They do not suffice to decide whether the solution can have a second regular center, i.e. perfect 
fluid solutions which occupy the whole ellipsoid. From equation ([5T|) we can conclude that the pressure is increasing 
near the first regular center if we assume A > 'iTTp{pc) + VZ-npc, however, we cannot control the further behavior of the 
solution and may obtain a singular solution where the pressure diverges, or a solution with a second center having a 
conical singularity. 

Let us furthermore discuss the consequences of having a regular center. Regularity of the solution at the center 
in particular fixes some of the coefficients in the power series expansion of the function R{r), see e.g. |46i [. which are 
given by 

i?(re) = 0, R'{r,) = l, i?"(r,)=0. (56) 

However, the part R"'{rc) is also determined by the initial conditions, as can be seen from the first field equation ([35]) 

, . , 1 — it 2ii , ^ 

87rp{pc) + A = lim — — lim — — . (57) 

1 — *r^ H r^rc H 



After applying the rule of L'Hopital and using the relations (|56p we arrive at 

8^p(Pe)+A = -3i?"'(re), (58) 

so that also the third derivative is fixed by the initial conditions. In case there exists a second regular center Tc^ we 
have 

i?(r,J=0, i?'(r,J = -l, R"irc,) = (59) 

so that R"'{rc2) is given by 

87rp(p,J+A = 3i?"'(r,,). (60) 

The conditions ([SG]) . ((58|) and also ([59|) . ((60|) must be satisfied independently by any solution admitting two regular 
centers. While one prescribes initial conditions at the first center, the regularity of the second center is by no means 
warranted since it actually depends on the solution of the function R together with the equation of state that also 
enter ([60]) . 

In the next section we will analyze the field equations numerically for given equations of state. It will turn out 
that none of the equations of state considered allows a second regular center. Therefore, all solutions that have 
an increasing pressure near the first center will have either a divergent pressure, or a second center with a conical 
singularity, and arc hence not of physical interest in general. 

IV. NUMERICAL CONSIDERATIONS 

In sections III Al and IIIDI we generalized two known solutions to include the cosmological constant. For large 
cosmological constant we found that new properties arise like the possibility of having a second regular center. In 
the previous analytic section we presented some arguments in favor of the existence of solutions with cosmological 
constant that may occupy more than 'half the three space. This result can also be read in the flowing way: For 
sufficiently large cosmological constants the pressure cannot vanish before the equator of the three space, so that some 
new physical properties may arise. However, it was also argued that in general a second regular center does not exist. 

It is the aim of the present section to solve the field equations (|35p - ([37)) numerically for a given equation of state 
p — p{p). In particular we are interested in solutions where the pressure vanishes after the group orbit's maximum 
of the three space and in solutions that possibly have two regular centers. We concentrate on these two classes of 
solutions, since the others are well known already. 

The most natural starting point for the numerical analysis are the polytropic equations of state 

p{p) = Kp^, p{p)^(l,)^\ (61) 

where K is some constant and n is the polytropic index. In the Newtonian case stellar models are finite if 1 < n < 5 
and do not have a finite radius for n > 5, where also A = is assumed. We slightly modify the polytropic equations 



of state, in order to allow a non- vanishing boundary density pb 
of state 
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p{p = 0). Hence, we consider the following equation 



P(P)-^ 



Pb, 



(62) 



where the boundary density is a new free parameter that we must specify. We chose K = \^ the polytropic index 
n = 3 and ph = 0.5. For two different cosmological constant ('small' and 'large') we obtain the two following solutions, 
Fig. Hi and Fig. Wp- 
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FIG. 4: Pressure function and radius for the polytropic equations of state, with K 
are pc = 1.0 {pc = 1.5) and A = 7r/2 (left), A = I.Stt (right). 



1, 71 = 3 and pb = 0.5. Initial conditions 



As expected, for small cosmological constant, see Fig. |4^, the pressure vanishes before the maximum of the group 
orbits. At the vanishing pressure surface one can join on the Schwarzschild-de Sitter metric C^ as the exterior 
spacetime, with the methods described in subsection IIIBI These solutions are represented by the Carter-Penrose 
diagram [1] discussed earlier. For large cosmological constant however, see Fig |4)d, the pressure vanishes after the 
maximum, in a region where the group orbits are decreasing. This would necessarily yield a global solution represented 
by the Carter-Penrose diagram[3l where the exterior spacetime contains the singularity. Since the numerical solutions 
vary smoothly in the cosmological constant, it is evident that a fine tuned cosmological constant can be chosen such 
that the pressure vanishes exactly at the maximum of the group orbits, in which case one has to join the Nariai metric 
as the exterior spacetime. 

Next, let us analyze the solutions for the stiff matter equation of state (the n ^ oo limit of the polytropic equation 
of state) 



p{p) =p + Pb, 



(63) 



which we, as before, supplemented by a boundary density term pb- For the stiff matter case we again take two different 
values of the cosmological constant. The results are similar to those discussed already, see Fig. [5^ and Fig.[S)D. 
Finally let us consider the Hagedorn equation of state 

P 



P{p) 



p expi 



1 



(64) 



where the free parameter p* is related to the boundary density by pb = p{p = 0) = p* /e. As in the previously discussed 
cases, we find that for 'small' cosmological constants the pressure vanishes before the maximum of the group orbits, 
whereas 'large' values of the cosmological constant allow the pressure to vanish after the maximum, see Fig. [6^ and 
Fig. B>- 

Apart from the incompressible perfect fluid case we could not find numerically any other configuration with a 
second regular center. This indicates that the existence of the second center requires a very specific choice of the 
fluid's equation of state and carefully chosen initial conditions. 



SUMMARY AND CONCLUSIONS 



We have used analytical and numerical techniques to analyze the static and spherically perfect fluid field equations 
of general relativity in the presence of a cosmological constant. A positive cosmological term can be viewed as an 
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FIG. 5: Pressure function and radius for the stiff matter equation of state, with p(, = 0.5. Initial conditions are pc = 1-0 
(pc = 1.5) and A = TT (left), A = S.Ttt (right). 
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external force having the effect of pulling matter apart. Hence, one can expect that the radial size of matter spheres 
is increased due to A. This naturally yields questions regarding the physical picture applicable to these solutions. It 
turns out that the effects of the cosmological constant lead to various different configurations, many of which have 
not been discussed previously. 

By using Gauss coordinates relative to the r = constant hypersurfaces, we analyzed geometrically the properties 
of the vanishing pressure surface that determines the boundary of the perfect fluid sphere. In the absence of the 
cosmological constant, going outwards, the area of the respective group orbits are always increasing close to the zero 
pressure surface. This situation changes drastically if A is allowed to be relatively large in comparison with the matter 
density. It is possible for the pressure to vanish exactly at the maximum of the group orbits or even vanish where the 
group orbits are decreasing. In the first case one has to join on the Nariai solution to get the metric C^ at boundary. 
In the latter case one matches the part of the Schwarzschild-de Sitter solution containing the black hole horizon and 
the singularity. This is in contrast to the the small A situation where the vacuum region contains the cosmological 
horizon. Lastly, one is lead to ask whether the matter can occupy the whole spacetime resulting in two regular centers 
corresponding to a fully generalized Einstein static universe where neither the energy density nor the pressure are 
constant. 

We showed that the Whittaker solution can have its vanishing pressure surface where the group orbits are decreasing, 
however, a second center is not possible. On the other hand, the Tolman IV solution does allow for a second regular 
center, a solution that might be named Tolman IV Einstein universe. By numerically integrating the field equations 
for physically motivated equations of state we showed that in general the pressure can vanish where the group orbits 
are decreasing and consequently at the maximum for sufficiently fine tuned initial conditions. We also obtained 
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analytical bounds that the cosmological constant has to satisfy to allow for such situations. However, we were not 
able to show, that in general, solutions with two regular centers for a given equation of state exist. This observation 
has analytical support since the conditions under which the solutions have two regular centers are very restrictive. 
Note however, that the special Tolman IV equation of state (pS)) admits solutions with a second regular center. 

Ever since the first exact matter solutions have been obtained, static and spherically symmetric perfect fluid 
spacetimcs have remained a subject of great interest. The presence of matter that effectively acts like a perfect fluid 
with unusual equations of state, such as p/ p ~ —1, drastically changes the geometry of known solutions. 
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APPENDIX A: FIELD EQUATIONS AND RIEMANN TENSOR 

The non-vanishing Riemann tensor components are 



-ly'^-iy", w^0^ 



1-R'- 



R^ r0 = i?'' r0 = TT, R 0t ~ R rjit = ^'^'"fT- (A-1) 



One can rewrite the field equations ([55)) - (|57|) to get 

R" . A 777 



^'^-^'P-U^- (A3) 



R ' 3 R^' 

277 

^3 3^ 



7.'^ + 7." = 47r(p + p)-5f-i (A4) 



and hence the Riemann tensor in terms of physical quantities 



R'''rt^~MP+P) + ^ + t (A5) 



i?3 3^ 





R'^04 = 


2m A 
i?3 + 3 ' 


.r0 

. r0 - 


= i?"^0 -- 


A 77T, 


net 

T- 0t ■ 


= R tj>t = 


A 777 



(A6) 

(A7) 
(A8) 
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